Looking for meson molecules in B decays 
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We discuss the possibility of observing a loosely bound molecular state in a B three-body hadronic 
decay. In particular we use the QCD sum rule approach to study a rj' — it molecular current. We 
consider an isovector-scalar I G J PC = 1 _ ++ molecular current and we use the two-point and three- 
point functions to study the mass and decay width of such state. We consider the contributions of 
condensates up to dimension six and we work at leading order in a 3 . We obtain a mass around 1.1 
GeV, consistent with a loosely bound state, and a 77' — % — > K + K~ decay width around 10 MeV. 
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One of the outstanding open questions in hadron physics is: are there meson-meson bound states? The same 
mechanism of meson exchange that binds the deuteron could also in principle bind two mesons. The interest in this 
subject was renewed by the discovery of the new charmonium states. Since their first appearance, some of them were 
considered to be meson molecules. These states have already been discussed in some reviews 

In this note we discuss how to look for molecules at the LHCb taking advantage of the unprecedented high statistics. 

We can look for meson molecules in three-body hadronic B decays. Since the phase space is large we can even try to 
use directly the Dalitz diagram, which extends up to large values of the variables S12 and S23. All the known normal 
quark-antiquark intermediate resonant states, leave an imprint in the Dalitz plot, which is directly related with the 
quantum numbers of the states and lead to the identification of the state. Examples are: a continuous straight line, 
in the case of scalar states, a line with a hole, in the case of vector states, or a line with two holes, in the case of 
tensor states. 




FIG. 1: Dalitz plot of a three-body B decay. The small drawings illustrate the different kinematical configurations. 



A sketch of a Dalitz plot for a three-body meson decay is shown in Fig. |H where for each invariant parameter S12 
or S23 it is shown the relative momentum of each one of the two particles 1, 2 or 2, 3. Let us consider the case that 
the particles 1 and 2 are pions coming from the p meson decay. Of course this decay should produce a line parallel 
to the S23 axis in the point s±2 = m 2 p . However, since the pions coming from the p meson decay must have one unit 
of angular momentum, they cannot go both to the same direction. Therefore, no pions could be seen in the region 
where the relative momentum between them is small. From Fig|T] one can see that this region is just in the middle 
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of the line parallel to the S23 axis. Therefore, a line characterizing a vector resonance state must have a hole in the 
middle, as mentioned above. Now imagine that the resonant state is a loosely bound molecular state of the particles 
1, 2. A loosely bound molecular state can only exist when the relative momentum between the two mesons in the 
molecule is small. In this case one has exactly the opposite situation than the one discussed before: there will be no 
signal in the Dalitz plot unless the two mesons in the molecular state go in the same diretion. Or, in other words, one 
expect a small line parallel to the S23 axis in the middle of the Dalitz plot, approximately in the region where there 
is a hole in the line characterizing a vector resonance. 

The final particles observed in the three-body B decays are pions and kaons. Therefore, to observe a molecular 
state in the Dalitz plot for a three-body B decay, this molecular state must decay into pions and/or kaons. Let us 
consider a rf — 7r loosely bound molecule with the quark content uuss. This resonant state, hereafter called R, is 
especially interesting because its mass, tur, should be approximately given by: 

m R ~ m v , + = 958 + 138 = 1096 GcV (1) 

and, therefore, it is quite visible in the B decay Dalitz plot. Since for S'-wave this molecule has I G J PC = 1~0 ++ , it 
can not decay into tt + tt~, but it will decay into K + K~ . In particular, there are already data for B~ — s- K + K~ K~ , 
B~ — > K + K~n~ , B~ —> 7r + 7r~7r~ and B~ —> tt + tt~K~ [4|. Only in the first two of these cases the decay could 
go through the resonant state i?, as illustrated in Figs. [5] and [31 For these cases, a small line with ^fsrz ~ 1.1 GeV 
parallel to the S23 axis should be seen in the Dalitz plot, in the region where the two particles 77' and it have a small 
relative momentum. This signal should be very different from all other established resonant states decaying into 
K + K~ , like the ao for instance, and should be only seen in the channels B~ — » K + K~K~ and B~ — > K + K~tt~ . Of 
course, the figure is very qualitative and it is not possible to say how large is the line segment around the indicated 
position. However, the observation of this structure in the Dalitz plot of the two mentioned B decays, and not in the 
others, would represent a strong evidence of the formation of this molecular state. The observation of a line that only 
appears in a certain piece of the S23 axis with a fixed value of S12, and only for the decays B~ — > K + K~ K~ and 
B~ — > K + K~n~ , could be interpreted as the existence of a weakly bound molecular state. 

II. THE 77' - 7T SCALAR MOLECULE 
A. Mass 

In a previous work [f| we have investigated the possiblity that the light scalar states could be interpreted as 
tctraquark states. Now we perform a complementary investigation, trying to understand a possible rj'ir meson molec- 
ular state in the QCD sum rule (QCDSR) framework [6j-|8j. 

The QCDSR approach is based on the correlator of hadronic currents. A generic two-point correlation function is 
given by 
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FIG. 3: The two relevant diagrams for the B — > K + K ir decay, through the resonance R. 



n(g) = i / d 4 x e^(0|T[j(x)^(0)]|0), 



(2) 



where the local current j(x) contains all the information about the hadron of interest, like quantum numbers, quarks 
content and so on. A molecular current can be constructed from the mesonic currents that describe the two mesons 
in the molecule. In the case of scalar rf — ir state a possible current is: 
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where i,j are color indices, u, d, s are the up down and strange quark fields respectively, and the mixing angle, 9, in 
the T]' current is 9 — 40° [a-fTlj. In this work we use 9 = 40°. 

In general, there is no one to one correspondence between the current and the state, since a molecular current can 
be rewritten in terms of a sum over tetraquark type currents through a Fierz transformation. However, as shown in 
, if the physical state is a molecular state, it would be better to choose a molecular type of current so that it has a 
large overlap with the physical state. In any case, it is very important to notice that since the current in Eq. @ is 
local, it does not represent an extended object, with two mesons separated in space, but rather a very compact object 
with two singlet quark- antiquark pairs. 

The coupling of the scalar resonance R, to the scalar current j, can be parametrized in terms of a parameter A as: 



(0\j\R) = A . 



(4) 



In the QCD side evaluation of the correlator function in Eq.([2]) we work at leading order and consider condensates 
up to dimension six. We deal with the strange quark as a light one and consider the diagrams up to order m s . We 
neglect the terms proportional to m u and m^. In the phenomenological side we consider the usual pole plus continuum 
contribution. Therefore, we introduce the continuum threshold parameter so [l2j . In the SU(2) limit the quarks u 
and d are degenerate and we consider the u-quark condensate equal to the d-quark condensate, wich we call (qq). 
After doing a Borel tranform in both sides of the calculation the sum rule is given by: 
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where 



n i i 

E =l-e- s »/ M T(-r- 



(6) 



fe=0 



which accounts for the continuum contribution. 

In the numerical analysis of the sum rules, the values used for the quark masses and condensates are (13 , flij : 
m a = 0.13 GeV, (qq) = '-(0.23) 3 GeV 3 , (ss) = 0.8(qq), (qga.Gq) = m%{qq) with m 2 , = 0.8 GeV 2 and (g 2 G 2 ) = 
0.88 GeV 4 . 

In Fig. ||]we show the OPE convergence of the sum rule in Eq. (JSJ) - From this figure we see that the convergence 
is reasonable for M 2 > 1.2 GeV 2 and very good for M 2 > 1.5 GeV 2 . However, as in the case of the light scalars 
14j, there is no pole dominance for these values of M 2 . This result could be interpreted in two different ways: i) it 
could indicate that this state does not exist, or ii) it could indicate that this state is not clearly separated from the 
continuum. The second interpretation can be applied to very broad states, as the light scalars a and k, since their 
widths are as large as the difference between their masses and the continuum threshold. In what follows we stick to 
the interpretation ii). 
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FIG. 4: The OPE convergence in the region 1.0 < M 2 < 2.0 GeV 2 for y/sE = 1.5 GeV. The dotted, dashed, long-dashed, 
dot-dashed, solid with circles and solid whit squares lines give the perturbative, quark condensate, gluon condensate, mixed 
condensate, four-quark condensate and total contributions repectively. 

In order to extract the mass ma without knowing the value of the constant A, we take the derivative of Eq. ([5]) 
with respect to 1/M 2 and divide the result by Eq. In Fig. [SJ we show the resonance mass as a function of M 2 for 
different values of y/so- We limit ourselves to the region M 2 > 1.2 GeV 2 where the curves are more stable and where 
the OPE convergence is better. Averaging the mass over all this region we find: 



m R = (1.15 ±0.10) GeV , 



(7) 



which is compatible with the experimental threshold in Eq. ([!]). Having the mass, we can also evaluate the value of 
the parameter A that gives the coupling between the state and the current. We obtain: 



A= (1.39 ±0.27) x 10" 3 GeV 5 



(8) 



B. Decay width 

In order to study the RK + K~ vertex associated with the R — > K + K~ decay, we consider the three-point function 



T^(p,p',«) = J d A x d A y e^'- x e*«(0\T{j£(x)j£-(y)j*(0)}\0), 
where p = p' + q, j is given in Eq. ([3]) and we use the axial currents for the kaons: 
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FIG. 5: The resonance mass as a function of the sum rule parameter (M 2 ) for different values of the continuum threshold: 
y'io = 1.4 GeV (dotted line), y/so = 1.5 GeV (solid line) and y/so = 1.6 GeV (dot-dashed line). 



To evaluate the phenomenological side we insert intermediate states for K + , K and R, and we use the definitions 
in Eqs. (0 and dTTJ bellow: 



<0|j& \K(p))=i Pll F K 



We obtain the following relation: 
T^ n {p,p',q) = 



F 2 X 



{M 2 R -p*){m 2 K -p> 2 ){m 2 K -q 2 ) 



9rkk v'ulv + higher resonances , 



where the coupling constant qrkk is defined by the matrix element: 

(K(p')K(q)\R(p))=g RKK . 



(11) 
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Here we follow refs. [jj UU and work at the kaon pole, as suggested in Hi for the nucleon-pion coupling constant. 
This method was also applied to the nucleon-kaon-hyperon coupling [ill Il7j. to the D* — D — n coupling [l8|, lilij 
and to the J/ip — it cross section [20l |. It consists in neglecting the kaon mass in the denominator of Eq. (fl"2l) in 
the term l/(m 2 K — q 2 ), and working at q 2 = 0. In the QCD side one singles out the leading terms in the operator 
product expansion of Eq.© that match the 1/q 2 term. Up to dimension six only the diagrams proportional to the 
quark condensate times m s and the four-quark condensate contribute. Making a single Borel transform to both 
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where sff = (1.0 ± 0.1) GeV 2 , is the continuum threshold for the kaon. 

As discussed in ref. [21| . the problem of doing a single Borel transformation, in a three-point function sum rule, 
is the fact that terms associated with the pole-continuum transitions are not suppressed. However, as shown in [2l| . 
the pole-continuum transition term has a different behavior, as a function of the Borel mass, as compared with the 
double pole contribution: it grows with M 2 . Therefore, the pole-continuum contribution can be taken into account 
through the introduction of a parameter A in the phenomenological side of the sum rule in Eq. (|14p , by making the 
substitution Qhk+k- ~^ 9rk+k- + am 2 [E [T3, EE 

Using Fk — 160 MeV, rajc = 490 MeV, wir — 1.15 GeV and the parameter A given by the sum rule in Eq. ([5]) we 
show, in Fig. [10 the QCDSR results for the vertex coupling constant, for different values of Sq and Sq in the interval 
given above. We see that, in the Borel range used for the two-point function, the QCDSR results do have a linear 
behaviour as a function of the Borel mass. Fitting the QCDSR results by a linear form: g^K+K- + AM 2 (which is 
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FIG. 6: The QCDSR result for the coupling constant g^K+K- as a function of the sum rule parameter M 2 for different 
values of so and So' (circles, triangles and squares for y^io = 1.5 GeV, s$ = 1.0 GeV 2 ; ^/io = 1.6 GeV, sff — 1.1 GeV 2 and 
y^so = 1.4 GeV, s^ = 0.9 GeV 2 respectively). The solid, dotted and dashed lines give the linear fit to the QCDSR results. 

also shown in Fig. |6]), the coupling can be obtained by extrapolating the fit to M 2 = 0. In the limits of the continuum 
thresholds mentioned above and taking into account the uncertainties in thr given in Eq. ((7]) we obtain: 

9rk+k- — (0.63 ± 0.06) GeV. (15) 

The decay width of R — > K + K~ is given in terms of the hadronic coupling gRK+K- as: 

T(R -+ K+K-) = ^ffW^K^I.O, (16) 

where X(m R , m 2 K , ra 2 K ) = m R + m^ K + — 2m 2 R m 2 K — 2m 2 R m 2 K — 2m 2 K m 2 K = rn 2 R {m R — Am 2 K ). Therefore, we get: 

T(R^K + K-) = (11.4 ±2.2) MeV. (17) 

Of course this is not the total width of the rfit molecule, since it can also decay into r\ — tt with a much bigger phase 
space. However, in the B decays discussed here, only the channel R — > K + K~ can be observed. 

The errors quoted above come directly from the uncertainty in the determination of the continuum threshold 
parameters, sq. According to our previous experience, they are the main source of uncertainty in the method. For a 
detailed analysis of the uncertainty associated to other parameters used in QCDSR we refer the reader to Refs. 
and 0. 

III. CONCLUSION 

We have proposed that a loosely bound molecular state should leave a particular signal in the Dalitz plot. A 
loosely bound molecular state, of the particles 1, 2, can only exist when the relative momentum between these two 
particles is small. Therefore, we expect to observe a short line parallel to the S23 axis in the middle of the Dalitz 
plot, approximately in the region where there is a hole in the line characterizing a vector resonance (see Fig. UJ). This 
signal is different from any signal charactering the normal quark-antiquark mesons, and could be used to identify the 
existence of loosely bound molecular states. 

In the case of three-body B decays, the final particles observed are pions and kaons. Therefore, to observe a 
molecular state in the Dalitz plot for a three-body B decay, this molecular state must decay into pions and/or kaons. 
We have considered a ?/ — tt molecular state. If this state exists as a loosely bound state, its mass should be close 
to the 77' — 7r threshold: ~ 1.1 GeV, quite visible in the B decay Dalitz plot. Since for a S'-wave this molecule has 
I G J PC = 1-0++ it can not decay into ir + ir , but it will decay into K + K . Therefore, the observation of a small line 
with y/Si2 ~ 1-1 GeV parallel to the S23 axis in the Dalitz plot for the B~ — > K + K~K~ and B~ —> K + K~n~ decays, 
with negative observation in the Dalitz plot for the B~ — > 7r + 7r~ K~ and B~ — > 7r + 7r~7r - decays would definitively 
indicate the existence of the rf — tt molecular state. 
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We have used QCD sum rules to study the mass and the decay width, of a 77' — tt molecular current, using two-point 
and three-point functions respectively. We have considered diagrams up to dimension six in both cases. We found 
a mass a slightly larger than the 77' — tt threshold, indicating the possiblity of a loosely bound molecular state. We 
obtained a small width for the 77' — tt — > K + K~ decay around 10 MeV. With these informations should be possible 
to experimentaly indentify this state in the B~ — > K + K~K~ and B~ — > K + K~tt~ Dalitz plots, if it exists. 

The method for the identification of resonances (or bound states) discussed here could be applied to other cases. 
A straightforward extension of our work could be done to the 77' — tt with quantum numbers J PC = 1 K This exotic 
state has been recently searched for by the CLEO [HJ and COMPASS (23[ collaborations. 
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